Measurement of the Generalized Polarizabilities of the Proton in Virtual Compton 
Scattering at Q 2 =0.92 and 1.76 GeV 2 : II. Dispersion Relation Analysis 
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Virtual Compton Scattering is studied at the Thomas Jefferson National Accelerator Facility in the 
energy domain below pion threshold and in the A(1232) resonance region. The data analysis is based 
on the Dispersion Relation (DR) approach. The electric and magnetic Generalized Polarizabilities 
(GPs) of the proton and the structure functions Pll — Ptt/e and Plt are determined at four- 
momentum transfer squared Q 2 = 0.92 and 1.76 GeV 2 . The DR analysis is consistent with the 
low-energy expansion analysis. The world data set indicates that neither the electric nor magnetic 
GP follows a simple dipole form. 

PACS numbers: 23.23.+x,56.65.Dy 



Virtual Compton Scattering (VCS) j*p — > jp has de- 
veloped in the last decade as a powerful tool to study 
the nucleon structure. At low CM-frame energy W, the 
VCS amplitude is parametrized as a function of the Gen- 
eralized Polarizabilities (GPs) of the proton which 
depend on the four-momentum transfer squared Q 2 of 
the virtual photon. These new observables have become 
the subject of experimental investigation via the study 
of photon electroproduction ep — > ep~/. In a companion 
Letter (referred to as I) we present an extraction of the 
structure functions Pt,l ~ Ptt/c and Plt from data be- 
low Nn threshold 13 , following the low-energy expansion 
(LEX) formalism Q. 

B. Pasquini et al. recently developed a formalism for 
VCS based on Dispersion Relations (DRs) ||. In this 
Letter we report a determination of the electric and mag- 
netic GPs of the proton (asiQ 2 ) and (3m{Q 2 ), respec- 
tively) from an analysis based on the DR model. The 
structure functions Pll — Ptt/c and Plt are also ex- 
tracted. The method uses the photon electroproduction 
cross section measured in the E93-050 experiment at 
Jefferson Lab (JLab). 

In the DR formalism, the VCS amplitude is predicted 
from the MAID parametrization [j| of pion electropro- 
duction, 7T° and er-meson f-channel exchange, plus two 
other Q 2 -dependent functions (subtraction constants). 
The latter are unconstrained phenomcnological contri- 
butions to otE{Q 2 ) and (3m(Q 2 )', therefore these GPs are 
not fixed by the DR model. They are written as: 

r exp ttN i 

a E (Q ) — a E (Q ) + - i - q2/a1 )2 (1) 

(same relation for (3m with parameter Ap) where a^ N 
{Plf) is the ttN dispersive contribution evaluated using 
the MAID analysis, and a e ^ p {(3 M ^ > ) is the experimental 
value at Q 2 — 0. The mass coefficients A Q and Kp are 
free parameters to be fitted experimentally. We note that 
the choice of a dipole form in Eq.^is not compulsory. A 
more fundamental property of the DR model is that, up 
to the Niti: threshold, it provides a rigorous treatment 
of the higher order terms in the VCS amplitude, beyond 
the lowest order GPs given by the LEX |l| . This feature 



allows the inclusion of data in the A(1232) resonance 
region in an extraction of GPs based on the DR approach. 

The JLab experiment uses the 4 GeV Continuous Elec- 
tron Beam Accelerator and the Hall A instrumenta- 
tionjl . More details can be found elsewhere 0, H, HI El 
EH Ip^ . The present study involves all the {ep — > epj) 
events with W < 1.28 GeV. The events are divided into 
three independent subsets listed in Table [I] For data 
sets I-a and II, by far most of the events lie below pion 
threshold; actually these two sets have also been used 
in the LEX analysis [2| considering only events such that 
W < (Mp + AU). Data set I-b covers mainly the A(1232) 
resonance region in W. The three data sets of Tableware 



TABLE I: Data sets for the DR analyses. 



data set 


Q 2 -range (GeV 2 ) 


VF-range 


I-a 


[0.85, 1.15] 


mostly < ttN threshold 


I-b 


[0.85, 1.15] 


mostly A resonance 


II 


[1.60, 2.10] 


mostly < ttN threshold 



the subject of three distinct DR analyses. The basic in- 
gredients for the cross-section determination are common 
to all our analyses of this experiment 0, H, 0, E3, 0, FlUI . 
They include a dedicated Monte-Carlo simulation [pj 
to obtain the acceptance, proper cuts to eliminate back- 
ground, the application of radiative corrections [T^j , the 
calibration of experimental offsets and the absolute nor- 
malization of the experiment. It is important to have 
a realistic shape for the sampling cross section in the 
Monte-Carlo in order to calculate an accurate solid an- 
gle. For this purpose, simulated events are sampled in 
the DR model cross section (d 5 <r DR ) which reproduces 
the enhancement of the A resonance. d 5 a DR depends on 
two free parameters A Q and A^ (cf. Eq. ^) which are it- 
eratively fitted by a x 2 minimization at the cross section 
level. When W increases, the acceptance is reduced to 
backward polar angles O-y-yCM (angle between the vir- 
tual and the final photons in the (<yp) CM frame). Cross- 
section data obtained in this angular region are shown 
in Figs, ^and [21 Kinematical conditions are defined on 
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Q 2 =1 GeV 2 , k lot , = 4.032 GeV , cos r , rC „= -0. 975 
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FIG. 1: (ep — > ep7) cross section for data set I-b in six in- 
tervals of the azimuthal angle <p (angle between lepton and 
hadron planes) as a function of W. By symmetry the statis- 
tics for cp= 180° to 360° are also included. Only statistical 
errors are shown. The solid curve is the prediction of the DR 
model for parameter values (A a , Ap) — (0.7, 0.6) GeV, while 
the dashed curve is for (0.5, 0.4) GeV and the dash-dotted 
curve for (0.9, 0.8) GeV, respectively. 



the plots, ki a b being the incoming beam energy (Fig.Hjl. 
q the virtual photon CM momentum and e the virtual 
photon polarization (Fig. |2J . Figure ^ clearly shows the 
A resonance excitation and the various curves indicate 
the sensitivity of the DR model to its free parameters. 
In contrast with Fig. ^ Fig- 121 shows how the Bcthc- 
Heitler+Born calculation, or the addition of a first-order 
GP effect as in the LEX approach 0, fails to reproduce 
the measured cross section above pion threshold. 




140 160 180/ -160 -140 140 160 180/ -160 -140 
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FIG. 2: (ep — » epy) cross section for data sets I-a (left) and II 
(right). The outgoing photon is emitted in the lepton plane 
(ip — 0° or 180°). The abscissa is the polar angle # 7 * 7 cm, 
with negative values corresponding to <p = 180° as in Ref. flo| . 
Only statistical errors are shown. The full curve is the DR 
model prediction using the values of (Aa, A/3) fitted to each 
of the two data sets (see Table ITU. The dashed (dash-dotted) 
curve is the BH+Born (plus a first-order GP effect) cross 
section. 

The results for Aq, and A^ are presented in Tabic ITT1 
Systematic errors are calculated from the same four un- 



certainties as in the LEX analysis |2( ■ The resulting error 
bars differ from one data set to another; this comes from 
the various data sets having a different phase space cov- 
erage, and a different sensitivity to both the physics and 
the sources of systematic errors. The reasonably good \ 2 
of the fits (1.3 to 1.5) indicates that the DR model works 
well in our kinematics and allows a reliable extraction of 
GPs, both below and above pion threshold. The compat- 
ibility between the three fitted values of A Q within errors 
(same for A^) suggests that the dipole form of Eq. ^is 
rather realistic, at least in the range Q 2 ~ 1-2 GeV 2 . 
We also note the close agreement between the obtained 
values of A a and A^. 

In the DR model these results directly translate into 
values for the electric and magnetic GPs, using Eq. ^ 
Table UTI gives the result of this evaluation of cye(Q 2 ) and 
Pm{Q 2 ) at Q 2 = 0.92 GeV 2 (data sets I-a and I-b) and 
Q 2 = 1.76 GeV 2 (data set II). These points are shown 
in Fig. |31 together with the point at Q 2 = llal and the 
points derived from the LEX analyses 0, Il5|. It must 
be noted that the latter do not directly yield the GPs, 
but the structure functions Pll — Ptt/^ and Plt, which 
are combinations of GPs. Therefore to obtain the "LEX 
points" of Fig. |21 we have extracted oir(Q 2 ) arid (3m{Q 2 ) 
from the measured structure functions Hid! according 
to the formulas 0: 

Pll - -Ptt = -^G V E a E {Q 2 ) + [spin-flip GPs] (2) 

£ tt QED 

2M I a 2 

Plt = P -\ j^G p E (3 M {Q 2 ) + [spin-flip GPs] (3) 

^QED V ^0 

where a QED is the fine-structure constant and G P E is the 
proton electric form factor. In this extraction the spin- 
flip GP terms are predicted by the DR model (so the 
result is model-dependent) and the parametrization of 
G P E is taken from Ref. 01 • The solid curve on Fig. |3 is 
the full DR calculation, split into its dispersive 7^ con- 
tribution (dashed curve) and the remaining asymptotic 
contribution (dash-dotted curve, dipole term of Eq.^) for 



TABLE II: Upper part: dipole mass parameters A a and A^ 
obtained by fitting the three data sets independently. Lower 
part: electric and magnetic GPs evaluated at Q 2 = 0.92 GeV 2 
(data sets I-a, I-b) and 1.76 GeV 2 (data set II). The first 
and second errors are statistical and total systematic errors, 
respectively. 



data set 


A a (GeV) 


Ap (GeV) 


I-a 


0.741 ± 0.040 ± 0.175 


0.788 ± 0.041 ± 0.114 


I-b 


0.702 ± 0.035 ± 0.037 


0.632 ± 0.036 ± 0.023 


II 


0.774 ± 0.050 ± 0.149 


0.698 ± 0.042 ± 0.077 


data set 


q b (Q 2 ) (10~ 4 fm 3 ) 


13m (Q 2 ) (10~ 4 fm 3 ) 


I-a 


1.02 ± 0.18 ± 0.77 


0.13 ± 0.15 ± 0.42 


I-b 


0.85 ± 0.15 ± 0.16 


0.66 ± 0.11 ± 0.07 


II 


0.52 ± 0.12 ± 0.35 


0.10 ± 0.07 ± 0.12 
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FIG. 3: Compilation of the data on electric (a) and magnetic 
(b) GPs. Data points are from Refs. |l6| (a), the LEX anal- 
yses of MAMI q3] (o) and JLab (o) and the present DR 
results (•). Some JLab points are shifted in abscissa for bet- 
ter visibility. The inner error bar is statistical; the outer one 
is the total error (quadratic sum of statistical and systematic 
errors) . The curves show the results of calculations in the DR 
model (see text). 



A«=0.70 GeV and A /3 =0.63 GeV, as fitted on the JLab 
data set I-b. The ttN intermediate states contribute very 
little to the electric polarizability at finite Q 2 (Fig. |3J-a) 
and they create a small dia-electric effect [a.^ < 0). 
The 7r7V contribution to the magnetic polarizability in 
Fig. [3jb is strongly paramagnetic, predominantly arising 
from the A(1232) resonance. In the DR formalism, this is 
cancelled by a strong diamagnetic term originating from 
tr-meson i-channel exchange, and parametrized by A^. 
The dotted curve is the full DR calculation evaluated for 
A Q = 1.79 GeV and A^=0.51 GeV, which reproduces the 
MAMI LEX data. The fact that there is no unique value 
of (A Q , Ap) agreeing with all data points suggests that 
the dipole form of Eq. ^ although working well in the 
range 1-2 GeV 2 as mentioned above, is not valid over 
the entire range of Q 2 . This further suggests that the 
global behavior of the GP a_g:(Q 2 ) does not follow a sim- 
ple dipole form. 

Finally, the results of our DR analyses can be expressed 



in terms of VCS structure functions. Using the GP val- 
ues of Table ILTland the DR model prediction for Ptt, the 
structure functions Pll~ Ptt / 'e and Plt are determined 
from Eqs.|21and|21 Their values are reported in Tabic ITTT1 
To facilitate the comparison with the LEX results, the 
present determination is made at similar kinematics in 
Q 2 and e (see Table 1111(1 . The agreement with the val- 
ues obtained in the LEX analysis |2( is very satisfactory. 
The results from the DR analysis of the (I-b) data set, 
obtained in the A region, yields smaller error bars due 
to an enhanced sensitivity to the GPs. 

In summary we have analyzed the process ep — > epj 
at JLab using a Dispersion Relation approach. The VCS 
structure functions obtained in this approach are in good 
agreement with the ones extracted using the low-energy 
expansion. We performed the first determination of GPs 
by analyzing data in the A(1232) resonance region. This 
opens up new possibilities to extract GPs from experi- 
ments, especially at higher Q 2 . 
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TABLE III: VCS structure functions obtained by the DR 
analyses of the three data sets. The first and second errors 
are statistical and total systematic errors, respectively. 



data 


Q 2 


e 


Pll - Ptt A 


Plt 


set 


(GeV 2 ) 




(GeV -2 ) 


(GeV" 2 ) 


I-a 


0.92 


0.95 


1.70 


-0.36 








± 0.21 ± 0.89 


± 0.10 ± 0.27 


I-b 


0.92 


0.95 


1.50 


-0.71 








± 0.18 ± 0.19 


± 0.07 ± 0.05 


II 


1.76 


0.88 


0.40 


-0.087 








± 0.05 ± 0.16 


± 0.019 ± 0.034 
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